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Abstract 

Superfluorescent decay in a weakly disordered slab is considered. 
By solving the Maxwell-Bloch equations we show, that the cooperation 
number increases with increasing of the mean free path. Superfluores- 
cent impulse is emitted at a small angle to the slab, isotropically in a 
plane direction. 

PACS.42.20-Propagation and transmission in inhomogeneous me- 
dia. 

PACS.42.20G-Scattering, diffraction and polarization. 
Introduction. 

During the spontaneous emission of uncorrelated exited atoms, their 
dipole moments interact each other via the emitted electromagnetic field. 
If the lifetime of exited state is long enough, the radiation can induce the 
correlation between the group of large number N c (so called cooperation 
number) atoms. As a result of correlation, N c exited atoms emit as one 
macroscopic dipole. Thus, the decay of the uncorrelated exited atoms, which 
starts with the low level of emitted radiation during the developing of macro- 
scopic dipole, is continued by the emission of this macroscopic dipole. At this 
last stage the duration of spontaneous decay decreases by the factor A^" 1 and 
the maximum of intensity of emitted impulse increases by the factor AT? com- 
pare that parameters of spontaneous decay of a single atom |L| . The total 
process is called superffuorescence @ . 

In extended system the cooperation number depends on the geometry 
of the sample. For a example, in a well studied pencil-like sample it is 
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N c ~ p\ 2 L in one mode approximation 0,0. Here p is density of the active 
atoms. L is the length of the sample, which diameter is much larger than 
the wavelength of radiation A. Maximum of impulse intensity is ~ N 2 x 
where N is total number of the exited atoms. The factor ^- ^> 1 can be 
interpreted as a number of macroscopic dipoles in the sample. 

Correlation between the exited atoms in this case develops due to the free 
propagation of the electromagnetic field along the pencil-like sample. During 
the propagation, the field can be absorbed and re-emitted. These processes 
put an upper limit on the geometrical size L, and therefore on N c 0, |J . 

The case, when the superfluorescence occurs in a disordered system, was 
considered in 0] . It was shown that in the disordered slab of the thickness 
L"^ I (lis the mean free path of the radiation due to scattering by the static 
fluctuations of dielectric constant) the cooperation number is N c ~ P^^r- 
Here factor -y- is the length of diffusive trajectory of the radiation. Scattering 
prevents the escape of the radiation from a system with the dimensions larger 
than the mean free path, therefore the cooperation number increases with the 
increasing of disorder up to the upper limit where absorption and re-emission 
starts. The maximum N c is of order of that in a pure case. 

The cooperation number is a random quantity in the disordered system. 
Approximations adopted in [|J in order to obtain the analytical solution 
allows calculation only of the average value of the cooperation number. 

Here we consider the superfluorescence of the weakly disordered slab with 
the thickness L less than the radiation mean free path. We show that the 
cooperation number in this system is N c ~ pl\ 2 . The correlation develops 
due to the propagation of the field along slab. Fluctuations of the dielectric 
constant scatter radiation out of volume, which contains the active atoms, 
thus decreasing the cooperation number. Superfluorescent impulse is emitted 
at small angle to the slab and isotropically in a plane direction. 

Maxwell-Bloch equations. 

We consider scalar version of problem. The coupling between the polar- 
ization density | {e luJt P (r; t) + e~ luJt P* (r; £)}, the population difference den- 
sity AN (f; t), and the field \ {e iu)t E (r; t) + e~ iujt E* (f; t)} can be described 
by the classical Maxwell-Bloch equations. In this approach amplified spon- 
taneous emission noise is neglected, which is a good approximation for the 
superfluorescence . P (r; t) and E (r; t) are slowly time- varying complex 
quantities, uo = ck is the atomic frequency. 

First two Maxwell-Bloch equations have the form 
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P(f;t) 



l \jJL 



AN (f; t) E (r; t) 



(1) 



n 



d 

dt 



AN (r; t) = {P* (f; t) P (r; t) - P (r; t) P* (r; *)} . 



(2) 



Here \i is the electric dipole moment. 

It is assumed that the population inversion relaxation time and the de- 
phasing time are larger than the delay time of the superfTuorescent impulse 
to, which will be defined lately in the equation (9). We also neglect inhomo- 
geneous broadening. 

The population difference AiV (r; t) and the polarization p ^' t ' in the equa- 
tions (|IJ) and (0) are the components of the local Bloch vector averaged over 
scales smaller than the wavelength of the radiation A. We choose the initial 
conditions P (r, t = 0) = and P(f,t = 0) P* (P, t = 0) = p \fi\ 2 5 (r — f ), 
which describe the initially uncorrelated atoms. The line means averaging 
over the initial state, r — (p, z) , p is in-plane coordinate. The population 
difference and the polarization are confined in a region of the slab \z\ < h. 

The field wave equation for the slow time- varying component E (r; t) has 
the form 



Here e(r) = 1 + 5e (r) is the dielectric function of the medium which 
contains active atoms. 8e (r) causes scattering. 
Calculation of cooperation number. 

We consider scattering by the small fraction of the particles of the size less 
than A embedded into the active medium. This is experimentally relevant 
case and for the theoretical treatment this case has an advantage that I is 
equal to the transport mean free path. 

Averaged over random positions of the scattering particles equation for 
the field has the form of the equation (H) with substitution of dumping ± for 



To consider superfluorescence we neglect time derivative in equation for 
the field at the very beginning. This usual approximation means that the 
time of the escape of radiation - from the system is smaller, than the time 
of exchange energy between field and atoms. 




(3) 



5e(r) 0. 
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Solving equation (JTJ) we obtain equation for the field, averaged over the 
disorder (E (r; t)) 



I.A + k 2 + j\(E (r; t)) = 4%k 2 P (f; £ = o) + 47rz ^ ^ | dtAN (r,t) (E (f;t)) 
^ > o 

(4) 

To find population difference from equation (2) we need to solve equation 
(4) inside the slab. Below we consider the case when I > kL 2 . This allows 
to neglect dependence of AiV (r, £) on coordinate r. We assume that at 
initial moment AiV (r, £ = 0) = p , where p is density of the active atoms. 
The quantity (E (r; £)) (P* (r; £)) weakly depends on coordinate when mean 
free path is larger than I > kL 2 . This allows to neglect the dependence of 
AiV (r, £) during the decay. 

From equation (||) it follows that the time dependence of field is deter- 
mined by the quantity 

2 r 

X = — ~ 2 dtAN(t) (5) 

pUT J 

where r = J- — = \l ^ Tr ?i is the characteristic time of energy ex- 

u y 2-Kpu\n\ V 3o;pA 0,; 

change between the field and the atomic system ||. r"^ = 8 ^^J is the 
radiative decay time of a single atom. 

Performing Laplace transformation of equation (|4]) over x an d Fourier 
transformation over in plane coordinate p we obtain solution inside the slab 

2nk 2 L ' 



(E (q; s)) = —r *, fcA1 / dz'P (q, z'; t = 0) , (6) 

2 

OO 

where (E (q; s)) = J d 2 pexp (iqp) J dxexp (-ax) (r; X (*))>, 



p = \Jk 2 — g 2 , Jm (p) > . 

The inverse Laplace transformation of @ gives an expression 



2ttA; 2 



2/ 



-exp 



Lk 2 X 



2 p+f 



dz'P (£ £ = 0) . (7) 
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Performing the inverse Fourier transformation and taking into account 
the initial conditions for the polarization density, we obtain 



(E(f; X ))(P*{r;x)) 



L 



ip |//| 



k 3 



8kl 2 x 



exp (2klx) 



Expression (||) gives asymptotic for kl\ > 1- 

Equation (0) and expression (§) allow to solve for the population differ- 
ence. Neglecting time-dependence of slow varying pre-factor in expression 
(8) we obtain for the population difference 



AN(f,t) = ptanh 



tp-t 
2t n ' 



where = and cooperation number is equal 



c 



6p\ 2 L 



(9) 



(10) 



» t n is calculated by matching (||) and 
21 and assuming that at the beginning of 





WTt, 



It is of order of number of active atoms in a tube with cross section A 
and length I. 

The delay time to ~ r^ln 
d^) with initial case when \ 

collective decay ( kl\ ~ 1) population difference does not change. 

Let us note that expression (^) coincides with the standard expression for 
the population difference in Markovian theory of superfluorescence. 

To obtain the expression for emitted radiation we must solve equation (f|) 
for \z\ > ^ and given time- dependent population difference @. 

Intensity of radiation from a unite volume of slab can be represented as 



I(f;t) = I (r) 



21 sin 6 
T 1 + | sin 9 



cosh ^- 
cosh 



(11) 



Here Iq (f) is the intensity of the ordinary fluorescence of a unite volume 
with exited atoms density p at distance r. It is isotropic function I (r) ~ r~ 2 . 

The second factor in (O) accounts for the mean free path of radiation in 
the medium, sin# = — . 

r 

The last factor in ([□]) describes the angular and time dependences of 
superfluorescent impulse. 



5 



Conclusion. 

We calculated the cooperation number and the intensity of cooperative 
decay in weakly disordered slab. Expressions (10) and (11) are the main 
results of the paper. 

It is shown that volume of the cooperating regions is proportional to 
volume of a tube with length, which is equal mean free path, and of cross- 
section dimension of order of the wavelength. The maximum of radiation 
is emitted in a small angle to the surface of slab isotropically in a plane 
direction. 

While calculation of (10) and(ll) was performed under the condition I » 
hi? it is reasonable to assume, that the scattering limits the cooperation 
number until I ~ L. 
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